Fluctuation theorem in driven nonthermal systems with quenched disorder 
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We demonstrate that the fluctuation theorem of Gallavotti and Cohen can be used to characterize 
the class of dynamics that arises in nonthermal systems of collectively interacting particles driven 
over random quenched disorder. By observing the frequency of entropy-destroying trajectories, we 
show that there are specific dynamical regimes near depinning in which this theorem holds. Hence 
the fiuctuation theorem can be used to characterize a significantly wider class of non-equilibrium 
systems than previously considered. We discuss how the fluctuation theorem could be tested in 
specific systems where noisy dynamics appear at the transition from a pinned to a moving phase 
such as in vortices in type-II superconductors, magnetic domain walls, and dislocation dynamics. 

PACS numbers: 05.70.Ln,05.40.Ca,74.25.Qt 



I 

o 



> 
o 

(N 

0^ 
O 

ON 

O 



X 



The fluctuation theorem (FT) of Gallavotti and Co- 
hen has been described as a generalization of the second 
law of thermodynamics in systems outside the thermody- 
namic limit P, 13; S| • It relates the frequency of entropy- 
destroying, also sometimes called second-law-violating, 
trajectories to entropy-creating trajectories and is suffi- 
ciently general to apply to systems far from equilibrium. 
It has been demonstrated to hold analytically for a class 
of time-reversible dynamical systems 111, and has been 
verified numerically in many others [J, 01- Wang et al. 
provided the first experimental verification of the rela- 
tion by observing the fluctuations of a dielectric particle 
pulled by an optical trap It has also been shown that 
the FT holds in some driven non-thermal systems such 
as granular materials [3| and a ball moving in a Sinai 
billiard potential Q. 

The FT has not previously been examined in the 
class of nonthermal nonequilibrium systems consisting 
of collectively interacting particles moving over a ran- 
dom background, where noisy dynamics can occur near 
pinned to moving transitions. Examples of this type of 
dynamics include the motion of magnetic domain walls 
[3| , depinning of electron crystals in solid-state materi- 
als [lOl.fll]|. plastic deformations in driven superconduct- 
ing vortex matter (l^ . [isl . [l3 . [l5|. and colloidal particles 
moving over quenched disorder |l6l. [l7j. Closely related 
to these systems are the dynamics of interacting dislo- 
cations under a strain Typically, there is a regime 
near the onset of motion where the particle trajectories 
are strongly disordered and 1// or crackling noise arises. 
In this work, we show that the FT can be used to charac- 
terize the dynamical behavior of a general model of this 
type, indicating that the FT could be applied to a much 
wider range of nonequilibrium systems than previously 
considered and may hold in the general class of systems 
exhibiting crackling noise. 



We speciflcally examine the formulation of the FT 
given in Ref. 0. One of the main predictions of the FT 
is that the probability density function (PDF) of the in- 



jected power p{Jt) obeys the following relation: 
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where J^is the injected power, r is the duration of the 
trajectory, and Sr is some constant. If St varies such 
that Pr — t/ Sr is constant for r 3> F^, where are the 
microscopic time scales of the system, then we say that 
Poo represents an "effective temperature." Wang et al. 
experimentally measured the quantity in Eq. 1 from the 
trajectories of a colloid that was driven through a thermal 
system. In the system we consider, there is no thermal 
bath; instead, the particle trajectories are generated in 
the presence of an external drive, a random quenched 
background, and interactions with other particles. 

We consider colloidal spheres confined to two dimen- 
sions and driven with an electric field in the presence 
of randomly distributed pinning sites. This particular 
model system has been shown to exhibit the same gen- 
eral dynamical features, including plastic fiow and mov- 
ing crystalline phases [19||, observed in other collectively 
interacting particle systems driven over random disorder 
such as vortices in type-II superconductors [l2,[l3|; thus 
we beheve the behavior in our system will be generic to 
other systems of this type. Additionally, experimental re- 
alizations of this system permit the direct measurement 
of the particle trajectories We simulate a system of 
Nc colloids with periodic boundary conditions in the x 
and y directions, and employ overdamped dynamics such 
that the equation of motion for a single colloid i is 



All quantities are rescaled 
damping constant rj is set 
arises from random Lang 



, + f^ + f; + fd (2) 

to dimensionless units, and the 
to unity. The thermal force 
evin kicks with the properties 
{{^) = and (f|,(i)f|,(<')) = 2r]kBTS{t - t')6^j. The col- 
loid interaction force fy is given by the following screened 



Coulomb repulsion: fy = 

Here Ac is an adjustable coefficient 



f ^)e-4-..f,^-. 
ri(^j^ is the position 
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of vortex Uj = jr, — rj\ and Vij = (r^ — rj)/rij. The 
quenched disorder introduces a force fp which is mod- 
eled by Np randomly placed attractive parabolic pin- 
ning sites of strength Ap and radius rp — 0.45, f* — 

Ylk=ii~^p'^ik/rp)Q{rp - nk)rzk, where 9 is the Heavi- 
side step function. The driving force id — fd^ is a con- 
stant unidirectional force applied equally to all colloids. 
We initialize the system using simulated annealing in or- 
der to eliminate undesirable transient effects due to re- 
laxation, and then apply the driving force. The equations 
of motion are then integrated by velocity Verlet method 
for 10^ — 10^ simulation time steps, depending on Nc. 
The time step dt = 0.002. We compute the longitudinal 
and transverse diffusivities Da with a = x,y hy fitting 
{[{r,{t + At) - r,{t)) ■ = 2Z?„At. 

The injected power computed for a single colloid i over 
a time period of length t is given by: 

rt+T 

Jr = {d-Ms)ds (3) 

where represents the instantaneous velocity of colloid 
i. A particle which moves opposite to the direction of 
the driving force makes a negative contribution to the 
entropy. We measure Jr for a series of individual par- 
ticles in a single run and combine this data to obtain 
p(Jt). We identify Jt for a variety of t ranging from 
a minimum of 10 simulation time steps to roughly one 
tenth the duration of the entire simulation. 

We first consider a system at fd — 0.1 with no pin- 
ning but with finite thermal fiuctuations T = 3.0 at a 
colloidal density of p = 0.5. Figure [DJa) shows p( Jt) 
for T = 0.42, 1.02, 1.62, 2.22, and 2.82. Over this 
range, p{Jt) is normally distributed with a slight right- 
ward skew due to the applied drive, and for increasing 
T the distribution sharpens. Equation [1] is certain to 
be followed since p( Jr) = C exp (- ((J) - Jr) /2o-2) and 
hence log {p{Jr) / p{- Jr)) = {-l/2a^)[{{J)-Jr)^-{{J) + 
Jr)'^] oc Jr. The validity of Eq. 1 for this system is illus- 
trated by the linear fits in Fig. [Ijb). 

We next repeat the procedure used to obtain Fig. 1 in a 
system with fd — 0.34, no thermal fiuctuations (T = 0), 
and which contains pinning sites with Ap = 0.5. For 
fd < 0.33 the particles are pinned and there is no non- 
transient motion, as illustrated in Fig. [SJa). Just above 
the depinning transition at fd = 0.34, the particle motion 
persists with time and the trajectories are highly disor- 
dered as shown in Fig. [2l[b). Approximately one third 
of the colloids are pinned at any given time; however, 
all of the particles take part in the motion. In Fig. [3l[a) 
we plot the strongly non-Gaussian p{Jr) that appear in 
the absence of thermalization for r = 0.02, 4.02, 8.02, 
12.02, and 16.02 at fd = 0.34. The t = 0.02 curve, 
most closely representative of the instantaneous distri- 
bution, peaks at Jt = and is skewed in the positive 
direction by the applied drive. Despite the strong non- 
Gaussianity of the PDF's, the ratio of entropy-producing 
to entropy-consuming trajectories is in agreement with 
the fiuctuation theorem of Eq. [1] as shown in Fig. [3l[b). 
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Figure 1: Demonstration of the FT for driven thermal par- 
ticles in the absence of pinning at p = 0.5, T = 3.0, and 
fd = 0.1. (a) Probability density function p(Jt) of injected 
power for all observed trajectories with r = 0.42, 1.02, 1.62, 
2.22, and 2.82 (from center bottom to center top), (b) Fit to 
Eq.[I]for r = 0.42, 1.02, 1.62, 2.22, and 2.82, from bottom to 
top. Linearity at various r indicates agreement with the FT. 
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Figure 2: Colloid positions (dots) and trajectories (lines) dur- 
ing 40000 simulation time steps at p = 0.5, Ap = 0.5, and (a) 
fa = 0.27; (b) fa = 0.34. 



The quality of the fit to Eq. [T] for fixed p and Ap de- 
pends on fd- It is known from earlier studies that systems 
with depinning transitions can exhibit a number of dif- 
ferent dynamical regimes as a function of external drive, 
including a completely pinned phase where there is no 
motion, a stable filamentary channel phase just at depin- 
ning where a small number of particles move in periodic 
orbits [12|, chaotic fiow at higher drives when the fila- 
ments change rapidly with time and a dynam- 
ically recrystallized phase at even higher drives where 
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Figure 3: Demonstration of FT in a nonthermal system with 
quenched disorder, (a) p{Jt) for all observed trajectories at 
p = 0.5 with Ap = 0.5 and fa = 0.34 at r = 0.02, 4.02, 
8.02, 12.02, and 16.02 (from upper right to lower right), (b) 
Fit to Eq. [l] showing agreement with the FT. Bottom to top: 
r = 0.02, 4.02, 8.02, 12.02, and 16.02. 



the particle paths are mostly ordered [1J|. To quantify 
the quality of the fits to Eq. 1 we calculate the Pearson 
product-moment correlation coefficient r [23|, which is 
a measure of the linearity of the relation between two 
variables. In Fig. 4(a) we plot the mean dissipation 
r]{v)/fd versus fd for the system in Fig. 2(b), along with 
the corresponding longitudinal and transverse diffusiv- 
ities Dx and Dy. In Fig. 4(b) we show the value of 
r for varied fd and for all r < Tc{fd), where Tc{fd) — 
sup {t I r{fd, r) > 0.5 V r < t}. Agreement with the FT, 
indicated by r « 1, holds over the largest range of r at 
fd ~ 0.34, near the depinning threshold and coinciding 
with peaks in both Dx and Dy. Here, the colloids flow 
in plastic fluctuating channels, as shown in Fig. 2(b). At 
lower drives fd < 0.33, and Dy are much smaller, 
the motion in the system is very unstable, and the par- 
ticles flow only through short-lived filaments, as shown 
in Fig. 2(a). We are unable to determine whether the 
FT fails to hold for fd < 0.33 since our measurement in 
this regime is dominated by rare events and our statistics 
remain poor over computationally accessible time scales. 
For higher drives fd ^ 0.4, the colloids begin to form 
an ordered crystal structure similar to that found in vor- 
tex systems at sufficiently high driving and both 
Dx and Dy drop. The FT continues to hold for small r 
at increasing fd, with the maximum value of r at which 
r > 0.5 decreasing with increasing fd- On short time 



Figure 4: Limits of regime in which FT is verified in a sys- 
tem with p = 0.5 and fp = 0.5. (a) Solid curve: mean 
dissipation ri{v)/fd vs fd, relating colloid displacements to 
applied drive. Upper crosses: longitudinal diffusivity 
vs fd- Lower crosses: transverse diffusivity Dy vs fd- (b) 
Pearson product-moment correlation coefficient r of the fit 
log(p( Jt)/p( — Jt)) = mJr + 6 as a function of r and fd- Val- 
ues closer to 1 indicate better agreement with the FT. The 
FT holds over the largest range of r in the fluctuating plastic 
flow regime near fd ~ 0.34 illustrated in Fig. 2(b). 



scales, the particles experience a "shaking temperature" 
Ts which decreases as oc 1/ fd Due to particle- 

particle interactions, on longer time scales the particles 
are effectively caged in a co-moving reference frame and 
no longer undergo long time diffusion. As a result, the 
FT fails to hold on longer time scales. 

As described previously, the FT allows the definition of 
an "effective temperature" Pt^oo when sufficient entropy- 
destroying trajectories of duration exceeding the micro- 
scopic time scales of the system can be sampled. This 
necessarily involves a balance of time scales since the 
second law of thermodynamics guarantees that p( Jt < 
0) = as T ^ cx). In Fig.O we plot (3r versus r showing 
the existence of an asymptotic effective temperature /3oo 
in a nonthermal system with quenched disorder. When 
we vary the initial configurations of the particle positions 
by changing the random simulation seed, we consistently 
find an asymptotic value of j3oo ~ 0.015 for r > 15. 
Equivalently, this indicates that the slope of our fits ob- 
tained as in Fig. [3Ub) scales such that t/St reaches a 
constant value at large t. 

We do not observe significant variation in Poo with fd; 
however, as noted previously, we can only define an ef- 
fective temperature for those values of fd where the FT 
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Figure 5: Effective temperature /3t in a nonthermal system 
with quenched disorder at p = 0.5, Ap = 0.5, and fd — 0.345. 
We overlay ten curves, each representing one realization with 
a unique random seed. The asymptotic value /3oo ~ 0.015. 
Inset: /32o = (/3i5<t<25) calculated over a range of p. 

holds over a wide range of t, which Umits us to drives 
near the depinning threshold where plastic flow occurs. 
To compare /3oo across different p, we perform our mea- 
surement at fd = 1-03/c for each p, where fc is the de- 
pinning force at that value of p. This places us within 
the plastic flow regime for every p considered here. The 
inset of Fig. [5] indicates that there is an apparently linear 
increase in (3oo with increasing particle density saturating 
at jO « 0.6. 



Experimentally testing the FT theorem for systems 
of collectively interacting particles in the presence of 
quenched disorder could be done in several ways. For 
superconducting vortices, the particle trajectories could 
be measured directly using various imaging techniques 
[2H . Recently it has been demonstrated that a single 
vortex can be dragged through a sample [IS], making it 
possible to perform a vortex experiment analogous to the 
colloid experiment of Wang et a/ Q . The most straight- 
forward measurement would be to monitor the voltage 
fluctuations at a constant drive to measure the power 
dissipation. A similar approach could be used to study 
conduction fluctuations in Wigner crystal systems 

We have shown that the fluctuation theorem of 
Gahavotti and Cohen can be apphed to a nonthermal 
system of collectively interacting particles driven over 
random disorder in the nonlinear regime above the depin- 
ning threshold where the particles flow plastically. This 
result indicates that the fluctuation theorem may be gen- 
eralized to a wide class of systems exhibiting fluctuating 
dynamics near a transition from pinned to moving, in- 
cluding magnetic domain wahs, vortices in type-II super- 
conductors, and sliding Wigner crystals. It would also be 
interesting to apply this approach to analyze other non- 
thermal systems that exhibit similar crackling noise such 
as dislocation dynamics. 
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